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The precise measurement of the cross section e"'"e~ — > 7r"'"7r~(7) from threshold to an energy of 3 GeV, using 
events with Initial State Radiation (ISR) collected with the BABAR detector, is presented. The ISR luminosity is 
determined from a study of the leptonic process e'^e~ — > jj.'^ fi~^{'y), and the method is tested by the comparison 
with the next-to-leading order (NLO) QED prediction. The leading-order hadronic contribution to the muon 
magnetic anomaly calculated using the BABAR tttt cross section measured from threshold to 1.8 GeV is (514.1 ± 
2.2(stat) ± 3.1(syst)) x 10^^". Other results on ISR multihadronic cross sections from BABAR are presented. 
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1. Introduction 

The measurements of the e~^e~ — hadrons 
cross-section are used to evaluate dispersion in- 
tegrals for calculations of the hadronic vacuum 
polarization (VP). In particular, the hadronic 
contribution to the muon magnetic moment 
anomaly (ajj"'*) requires data in the low mass re- 
gion, dominated by the process e~^e~ — > 7r"''7r~(7) 
which provides 73% of the contribution. The 
dominant uncertainty comes also from the tttt 
channel. 

The systematic precision of the (recent) pre- 
viously measured tttt cross sections is of 0.8% 
for CMD2 and 1.5% for SND [3], the two 
measurements being in good agreement. The 
first ISR measurement done by KLOE [J] had a 
quoted systematic precision of 1.3%, and some 
deviation in shape was observed when compar- 
ing to the Novosibirsk data. For the reanalysed 
KLOE data IS] a systematic uncertainty of 0.9% is 
quoted, and the agreement with the Novosibirsk 
data is improved. 

The comparison of the theoretical and mea- 
sured [1] values of shows a discrepancy of 
about 3(7 when previous e+e~ data |2|3|5j are 
used, possibly hinting at new physics. An ap- 
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proach using r decay data corrected for isospin- 
breaking, leads to a smaller difference |Xi- 

The BABAR 27r results presented in these pro- 
ceedings were published in [8]. Their achieved 
goal was to obtain a measurement of the contri- 
bution of the 2tt channel to [a^"''^) with a precision 
better than 1%, implying a control of systematic 
uncertainties at the 10"'^ level. 



2. The BABAR ISR tttt analysis 

The results on the tttt cross section presented 
here are obtained with the ISR method [6] us- 
ing e^e~ annihilation events collected with the 
BABAR detector, at a center-of-mass energy ^/s 
near 10.58 GeV. We consider events e^e" — 
X^iSFt., where X includes any final state, and 
the ISR photon is emitted by the e+ or e~ . The 
e+e^ — > 'm^{'^FSR) cross section is obtained as 
a function of \/s', which is the invariant mass 
of the final state. The advantage of the ISR 
method is that all the mass spectrum is covered 
at once (from threshold to 3 GeV in BABAR) 
with the same detector conditions and analysis. 

In the BABAR analysis the t^~^tt~ is ail fsr) 
and ^J'~ 1 1 s r{i F s r) spectra are measured. 
This is the first NLO measurement, an eventual 
additional radiation being taken into account in 
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the analysis, instead of being corrected a posteri- 
ori (as done by previous experiments). The muon 
spectrum is compared with the NLO QED predic- 
tion, this important cross check of the analysis 
being called the QED test. The VT' spectrum of 
e+e~ — Xj events is related to the cross section 
for the process e'^e~ X through 
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where ex-y is the detection efficiency (acceptance) 
determined by simulation with corrections ob- 
tained from data, and ct^ is the bare cross sec- 
tion (excluding VP). The muon spectrum is used 
to derive the effective ISR luminosity Lffj^. We 
correct for the leading order FSR contribution 
for muons (smaller than 1%, below 1 GeV), 
while additional FSR photons are measured. The 
TTTT^-jpsfi) cross section is obtained from the ratio 
of the 27r spectrum and Ljgj^, where the e+e" lu- 
minosity, additional ISR effects, vacuum polariza- 
tion and ISR photon efficiency cancel, hence the 
strong reduction of the systematic uncertainty. 

This analysis is based on 232 fb~^ of data 
recorded with the BABAR detector [9, at the 
PEP-II asymmetric-energy e~^e~ storage rings. 
Charged-particle tracks are measured with a five- 
layer double-sided silicon vertex tracker (SVT) 
together with a 40-layer drift chamber (DCH) in- 
side a 1.5 T superconducting solenoid magnet. 
The energy and direction of photons are mea- 
sured in the CsI(Tl) electromagnetic calorime- 
ter (EMC). Charged-particle identification (PID) 
uses ionization loss dE/dx in the SVT and 
DCH, the Cherenkov radiation detected in a ring- 
imaging device (DIRC), and the shower deposit 
in the EMC (Ecai) and in the instrumented flux 
return (IFR) of the magnet. 

Two-body ISR events are selected by requiring 
a photon with E* > 3 GeV and laboratory po- 
lar angle in the range 0.35 — 2.4 rad, and exactly 
two tracks of opposite charge, each with momen- 
tum p > 1 GeV/c and within the angular range 
0.40 — 2.45 rad. Events with one single charged 
track are also recorded and exploited for efficiency 
measurements. 

Signal and background ISR processes are sim- 
ulated with Monte Carlo (MC) event genera- 



tors based on Ref. [10]. Additional ISR pho- 
tons are generated with the structure function 
method [11], and additional FSR photons with 
PHOTOS [12]. The response of the BABAR de- 
tector is simulated with GEANT4 [14] . 

Background events from e+e^ ~^ 99 (<7 = 
u,d,s,c) are generated with JETSET [13]. They 
are due to low-multiplicity events in which an en- 
ergetic 7 originating from a tt^ is mistaken as the 
ISR photon candidate. To normalize this rate 
from JETSET, the tt'^ yield (obtained by pairing 
the ISR photon with other photons in the event) 
is compared in data and MC. Multi-hadronic 
ISR backgrounds are dominated by e~^e~ — > 
7r+7r"7r°7 and e+e~ 7r+7r^27r°7 contributions. 
An approach similar to that for qq is followed to 
calibrate the background level from the Stt ISR 
process, using u and cj) signals. The MC esti- 
mate for the 27r27r°7 process is used and assigned 
a 10% systematic uncertainty. Background con- 
tributions to the /i/i spectrum are found to be 
negligible. 

Acceptance and mass-dependent efficiencies for 
trigger, reconstruction, PID, and event selection 
are computed using the simulation. The ratios of 
data and MC efficiencies have been determined 
from specific studies, as described below, and are 
applied as mass-dependent corrections to the MC 
efficiency. They amount to at most a few percent 
and are known to a few permil level or better. 

Tracking and PID efficiencies are determined 
taking advantage of pair production. For tracking 
studies, two-prong ISR candidates are selected on 
the basis of the ISR photon and one track. A kine- 
matic fit yields the expected parameters of the 
second track. The unbiased sample of candidate 
second tracks is used to measure track reconstruc- 
tion efficiency. The study of 2-particle overlap in 
the detector required a large effort to reach the 
per mil accuracies. 

Each event is subjected to two kinematic fits 
to the e+e^ — > X'y hypothesis, where X allows 
for possible additional radiation. Both fits use 
the ISR photon direction and the parameters and 
covariance matrix of each charged-particle track. 
The two-constraint (2C) 'ISR' fit allows an un- 
detected photon collinear with the collision axis, 
while the 3C 'FSR' fit is performed only when 
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an additional photon is detected. Most events 
have small values for both fits; an event with 
only a small xjsR (Xfsr) indicates the pres- 
ence of additional ISR (FSR) radiation. Events 
where both fits have large values result from 
track or ISR photon resolution effects, the pres- 
ence of additional radiated photons, or multi- 
hadronic background. To accommodate the ex- 
pected background levels, different criteria in the 
(XisR'Xfsr) plane are applied, i.e. a loose 2D 
cut in the central p region and a tighter cut for 
the p tails. The loose cut is also used in the /i/17 
analysis. The tttt {pp) mass is calculated from 
the best TSR' or 'FSR' fit. 

The computed acceptance and the selection 
efficiency depend on the description of radiative 
effects in the generator. The difference of the FSR 
rate between data and MC is measured, resulting 
in a small correction for the cross section. More 
significant differences are found between data and 
the generator for additional ISR photons, since 
the latter uses a coUinear approximation and an 
energy cut-off for very hard photons. Induced 
kinematical effects have been studied using the 
NLO PHOKHARA generator [TS] at four-vector 
level with fast simulation. Differences occuring 
in acceptance yield corrections to the QED test. 
In contrast, since radiation from the initial state 
is common to the pion and muon channels, the 
7r7r(7) cross section, obtained from the tttt/ pp ra- 
tio, is affected and corrected only at a few per- 
mil level. The selection efficiency determined 
from muon data applies to pions, after correction 
for the effect of secondary interactions and the 
tt/p difference for additional FSR. The measured 
7r7r(7) cross section is to a large extent insensitive 
to the description of NLO effects in the generator. 

3. The QED test 

The QED test involves two additional factors, 
both of which cancel in the tttt /pp ratio: L^e and 
the ISR photon efficiency, which is measured us- 
ing a ppj sample selected only on the basis of the 
two muon tracks. The QED test is expressed as 
the ratio of data to the simulated spectrum, after 
the latter is corrected using data for all known 
detector and reconstruction differences. The gen- 
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Figure 1. (a) The ratio of the measured cross section 
for e"'"e~ — >■ p^ P~'y{'y) to the NLO QED prediction. 
The band represents Eq. ([2]). (b) The measured cross 
section for e^e~ — >■ 7r-'-7r~(7) from 0.3 to 3 GeV. (c) 
Enlarged view of the p region in energy intervals of 
2 MeV. The plotted errors are from the sum of the 
diagonal elements of the statistical and systematic co- 
variance matrices. 



erator is also corrected for its known NLO de- 
ficiencies using the comparison to PHOKHARA. 
The ratio is consistent with unity from threshold 
to 3 GeV/c^ (Fig. [1] (a)). A fit to a constant 
value yields (x^/fidf = 55.4/54; ndf=number of 
degrees of freedom) 



data 
^^7(7) 



NLO QED 



-1 = (40±20±55±94)xl0-'',(2) 



where the errors are statistical, systematic from 
this analysis, and systematic from Lee (mea- 
sured using Bhabha scattering events), respec- 
tively. The QED test is thus satisfied within an 
overall accuracy of 1.1%. 

4. Angular distribution in the tttt center- 
of-mass 

The pion angular distribution in the tttt center- 
of-mass with respect to the ISR photon direction 
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Table 1 

Relative systematic uncertainties (in 10^^) on the 
e+e~ — >■ 7r+7r~(7) cross section by \/s' intervals (in 
GeV) up to 1.2 GeV. The statistical part of the effi- 
ciency uncertainties is included in the total statistical 
uncertainty in each interval. 
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Figure 2. The angular pion distribution in the 
Tin system with respect to the ISR photon direc- 
tion as function of | cos 9^ \ for background-subtracted 
7r7r7(7) data (points) in the p central region (0.5 < 
m^^ < 1 GeV /(?). The blue histogram is the shape 
obtained in the simulation, normalized to the data. 



in that frame, is model-independent. The cos0* 
distribution behaves as sin^ 0* as a consequence 
of the P-wave between the two pions, but it is 
strongly distorted at | cos 0* | values near one due 
to the p> 1 GeV/c cut on the tracks. 

The I cos 6** | distributions for background- 
subtracted data and MC are compared in Fig. [2] 
for the 0.5 — 1 GeV/c^ mass range: they agree 
with each other within the statistical uncertain- 
ties, as expected for a pure pion sample. 

5. The TTTT cross section 

To correct for resolution and FSR effects, an 
unfolding of the background-subtracted m^jr dis- 
tribution, corrected for data/MC efficiency dif- 
ferences, is performed. A mass-transfer matrix, 
created using simulation, provides the probabil- 
ity that an event generated in a ^fs' interval i 
is reconstructed in a m^n^^r interval j. The ma- 



trix is corrected to account for the larger fraction 
of events with bad values (and consequently 
poorer mass resolution) in data compared to MC 
because of the approximate simulation of addi- 
tional ISR. The performance and robustness of 
the unfolding method [H] have been assessed us- 
ing test models. 

The results for the e+e^ — > tt+tt^ (7) bare cross 
section including FSR, cr°^^^^('\/s'), are given in 
Fig. [T] (b). Prominent features are the dominant 
p resonance, the abrupt drop at 0.78 GeV due to 
p—Lo interference, a clear dip at 1.6 GeV resulting 
from higher p state interference, and some addi- 
tional structure near 2.2 GeV. Systematic uncer- 
tainties are reported in Table [T] for 0.3 < ^/s' < 
1.2 GeV. Although larger outside this range, they 
do not exceed statistical errors over the full spec- 
trum for the chosen energy intervals. In partic- 
ular, a systematic uncertainty of only 0.5% has 
been achieved in the central p region. 

The BABAR data were compared to other ex- 
periments, exploiting a VDM fit of the pion form 
factor [IT]- This fit describes well the BABAR 
data in the region of interest for the compari- 
son. There is a relatively good agreement when 
comparing to the Novosibirsk data |2|3) in the 
p mass region, while a slope is observed when 
comparing to the KLOE '08 data 5 . A flat- 
ter shape is observed when comparing to the re- 



5 



cent KLOE [18 data, obtained by the analysis 
of events with a detected, large angle ISR pho- 
ton. A good agreement is observed when com- 
paring to the Novosibirsk and KLOE data, in 
the low mass region, below 0.5 GeV/c^. There 
is a good agreement between the BABAR data 
and the most recent r data (corrected for isospin- 
breaking effects) from Belle [19], while some sys- 
tematic effects are observed when comparing to 
ALEPH pni and CLEO [H]. 

6. The TTTT contribution to 

The lowest-order contribution of the 7r7r(7) in- 
termediate state to the muon magnetic anomaly 
is given by 

oo 

where K{s') is a known kernel [22]. The inte- 
gration uses the measured cross section and the 
errors are computed using the full statistical and 
systematic covariance matrices. The systematic 
uncertainties for each source are taken to be fully 
correlated over the full mass region. The inte- 
grated result from threshold to 1.8 GeV is 

^7r,r(7),LO ^ (5^4 ^ 2.2 ± 3.1) X 10"^" , (4) 

where the errors are statistical and systematic. 
This value is larger than that from a combination 
of previous e+e" data [7] (503.5 ± 3.5), but is in 
good agreement with the updated value from r 
decay [7] (515.2±3.4). The deviation between the 
BNL measurement [1] and the theoretical predic- 
tion is reduced to 2.4 standard deviations, when 
using the tt+tt" data from BABAR only. 

7. Other ISR measurements from BABAR 

Many other ISR measurements 23|24l25|26l27) 

have been published by BABAR in the past (see 
Fig- IS]) • They have an improved precision compar- 
ing to previous measurements. Still more chan- 
nels are under analysis: K'^ K~ , KKirn including 
K'^'s, and 7r+7r^27r°. The last one is especially 
important in order to clarify the existing discrep- 
ancy between the e^e~ and r-based contributions 
to from this channel [5S]. 
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Figure 3. Cross sections obtained by BABAR using 
the ISR method. The error bars include statistical 
and systematic uncertainties added in quadrature. 



8. Conclusions and perspectives 

BABAR has analyzed the tt'^tt^ and fJ,^ ^J,~ ISR 
processes in a consistent way, on the full mass 
range of interest (0.3 - 3 GeV/c^). In addition, 
the absolute /i^/i^ cross section has been com- 
pared to the NLO QED prediction, the two be- 
ing in agreement within 1.1%. The e+e^ — > 
7r+7r^(7) cross section, obtained through the ra- 
tio of the 7r+7r~ and fJ-'^ fJ.~ spectra is rather in- 
sensitive to the detailed description of radiation 
in MC. A strong point of the present analysis, 
comparing to previous ISR studies, comes from 
the fact that several uncertainties cancel in this 
ratio, allowing us to achieve our precision goal: 
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the systematic uncertainty in the central p region 
(0.6 - 0.9 GeV/c2) is only 0.5%. 

The contribution to from the BABAR tt+tt^ 
spectrum, in the range 0.28 — 1.8 GeV, is (514.1 ± 
2.2 ± 3.1) X 10"^". This result has a precision of 

0. 7., comparable to the combined previous mea- 
surements. The contribution from multi-hadronic 
channels will continue to be updated, with more 
results forthcoming from BABAR. 

In the comparison between the BABAR 
7r+7r^(7) cross section and the data from other 
experiments, there is a fair agreement with CMD2 
and SND, while the agreement is poor when com- 
paring with KLOE. The first priority should be to 
clarify the BABAR/KLOE discrepancy, the most 
important effect on being due to the differ- 
ence on the p peak. The origin of the slope in 
this comparison is also to be understood. The 
slope was very pronounced when comparing with 
the 2004 KLOE results, and it is reduced with 
the 2008 and 2010 data. The same slope is also 
observed in the comparison of the KLOE and r 
data, while BABAR is in good agreement with 
the T results. Further checks of the KLOE results 
are possible. Actually, as the method is based on 
MC simulation for ISR and additional ISR/FSR 
probabilities, an important cross-check should be 
provided by the comparison of the /x/i spectrum 
with the QED prediction. 
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